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Abstract 

We apply the MC@NLO formalism to the production of heavy-quark pairs in 
pointlike photon-hadron collisions. By combining this result with its analogue 
relevant to hadron-hadron collisions, we obtain NLO predictions matched to 
parton showers for the photoproduction of QQ pairs. We compare MC@NLO 
results to the measurements of c- and 6-flavoured hadron observables performed 
by the HI and ZEUS collaborations at HERA 
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1. Introduction 

The production of pairs of quarks with mass much larger than the typical 
hadronic scale, thus called heavy, has distinctive features that renders it an in- 
teresting case study. From the theoretical viewpoint, the quark mass cuts the 
collinear singularities off, and open-quark cross sections, in which no fragmen- 
tation functions or jet-finding algorithms are used, are well defined in perturba- 
tion theory. On the experimental side, heavy-quark production can be exploited 
thanks to its peculiar signatures, and generally large rates. The definition of 
what is heavy is to a certain extent ambiguous, and may depend on the kine- 
matic region one is interested into probing (since the behaviour of a particle 
with mass m can be fairly different according to whether m > or m < Pt)- 
While there is no doubt that the top is heavy, the heaviness of the bottom quark 
may be debatable, and that of the charm quark certainly is. It remains true, 
however, that open charm and bottom cross sections can be computed at fixed 
order in QCD; the interesting question is therefore how predictions compare to 
data. 
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This question has indeed received a lot of attention in the past twenty years, 
in part because of the important role played by heavy quarks in many discov- 
ery channels of BSM physics. Apart from the mainstream activity at hadron 
machines, the ep collider HERA has also been collecting charm and bottom 
data for over a decade now, in both the DIS and the photoproduction regimes 
(see e.g. ref [l| for a review). With its relatively small cm. energy of about 
300 GeV, heavy-quark physics at HERA very rarely involves large transverse 
momenta, and therefore offers a good testing ground for purely-perturbative 
QCD predictions. It is now customary to compare data to results accurate to 
NLO (C(a em a|)) [l|; fully exclusive computations, in both DIS @ and photo- 
production 0], have been available for some time now. 

Comparisons between data and parton-level, fixed-order predictions are not 
entirely satisfactory. They may work well for fully inclusive observables dom- 
inated by hard scales, but unfortunately experimental measurements typically 
involve hadron-level, smallish-scale quantities which cannot be well described by 
such a simple theoretical framework, and corrections (which introduce biases) 
need be applied to theoretical results, data, or both, in this way blurring the 
picture. Parton Shower Monte Carlos (PSMCs) offer a viable alternative, with 
their fully realistic final states, but lack the accuracy of higher-order perturba- 
tive computations. The solution to the problem of matching the two approaches 
has been extensively studied in the past ten years, and is now fairly well estab- 
lished (see e.g. ref. [4J for a pedagogical introduction). 

The purpose of this paper is that of applying the MC@NLO matching for- 
malism [5[ to the case of heavy-quark pair photoproduction, and of comparing 
the (hadron-level) results obtained in this way with selected bottom and charm 
measurements performed by the HI and ZEUS collaborations at HERA. This 
paper is organized as follows: in sect.[5]we give the briefest introduction to the 
photoproduction jargon; in sect.[3]we summarize a few technical information on 
MC@NLO, specific to photon-initiated processes. Section [4] presents the com- 
parisons between theoretical predictions and data; finally, in sect. [5] we report 
our conclusions. 

2. Photoproduction: generalities 

The high-energy photons that initiate photoproduction processes are typi- 
cally obtained by quasi-collinear bremsstrahhmg off an electron beam. We write 
the factorization theorem relevant to photoproduction as follows: 

dcr eH = ^1 dxidx 2 fif ) {xi)f { b H) (x 2 )do lb {xi,X2) 
b 

+ X] / dx i dx 2fa e H x l)fb H) ( X 2) d(J ab{xi,X2). (1) 
ab 

We have denoted by f-f 1 the Weizsiicker- Williams function 0, 0] , which gives a 
good approximation of the energy spectrum of photons radiated by the electrons 
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The functions /i e ' , the electron PDFs, are the analogues of the hadron PDFs 



(for the computations of this paper, we have used the form presented in ref. [8|). 
The functions /i e ' , the ele 
and are defined as follows: 

fi e \x) = I dydzS(x - yz)&\y)tt\z) , (2) 



with /a 7 ' the relevant photon PDFs. The two terms on the r.h.s. of eq. ([1} arc 
usually called the pointlike and the hadronic photon components respectively. 
As the name suggests, the hadronic component is computed in exactly the same 
way as any ordinary hadron-hadron cross section. MC@NLO is no exception 
to this rule, and we have therefore used the calculation of ref. [9], with trivial 
modifications in the corresponding computer code that allow the use of electron 
PDFs. In this work, we have performed the computation of the pointlike con- 
tribution to QQ production in MCQNLO. The underlying partonic processes 
are: 

1 + 9 — > Q + 0, (3) 
1 + 9 — > Q + Q + 9, (4) 
1 + q — > Q + Q + q, (5) 

where eq. Q receives contributions from the Born (0(a cm a s )) and one- loop 
(0(a cm a 2 s )) matrix elements, and eqs. (J4j) — (J5J) receive contributions from the 
real-emission matrix elements (0(a em ctg)), and from the corresponding coun- 
terterms as defined by the subtraction formalism adopted for the pure-NLO 
computation Q. The matrix elements have been taken from ref. Q, which we 
have subsequently matched to fortran HERWIG 1(1 U, 12 1, according to the 



MC@NLO formalism as described in the following section. 



3. MCONLO for QQ photoproduction 

The MC@NLO formalism has been introduced in ref. and applied since 
then to a fairly large number of hadroproduction processes. The relevant tech- 
nical details can be easily found in the literature, and we shall therefore refrain 
from giving them again here. We limit ourselves to recall that in the context of 
MC@NLO the matching of an NLO computation with a PSMC requires one to 
modify the short-distance cross sections that enter the former, with the inclu- 
sion of the so-called Monte Carlo (MC) subtraction terms, that are responsible 
for removing any double counting at the NLO. In turn, the MC subtraction 
terms have a factorized structure, being essentially constructed with Born-level 
cross sections and with process-independent branching kernels, which can be 
computed once and for all once the PSMC is chosen that will be used in the 
shower phase. 

In order to determine the right combination of the Born matrix elements and 
branching kernels that enter the MC subtraction terms for a given production 
process, one formally expands the all-order PSMC cross sections to NLO. The 
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way in which the results so obtained are manipulated to construct the MC sub- 
traction terms used in computer programmes is straightforward, and has been 
described in several papers. Here, we only give the results for the perturbative 
expansion mentioned above. Given the similarity between heavy-quark photo- 
and hadroproduction, we adopt the same notation as in ref. 0] , where the latter 
process was studied. The PSMC cross sections read: 



da 



(L,Z) 
eh 



(6) 



where the first sum in eq. ([6]) runs over real-emission parton processes, eqs. Q 
and |5]). The index L runs over the emitting legs and it assumes the values +, 
— , Q, and Q. The index I runs over the colour structures, and it can take the 
values t and u (clearly, there is no s-channel colour connection in the case of 
photoproduction). Following ref. [§], we obtain 



j (+,0 


MC 


z\> 




dxu 

MC 


dx 2i , 


(7) 


j (-.0 

d<J eb 


MC 




da^ 


dxu 

MC 


dx 2i , 


(8) 






= rt e \x lf )fl H \ X2f )d*^ 


dx\f 


dx 2f , 




(9) 


MC 


MC 








dog* 


MC 


= f^{x lf )fi H \x 2f )da^ 


dx\f 

MC 


dx 2 f . 




(10) 



These equations are identical to eqs. (5.2)-(5.5) of ref. [9j, apart from the obvious 
notational changes due to the incoming electron in place of a hadron. The short 
distance cross sections in eqs. (f7])-(fT0|) are: 
• 75 initial state (I — t, u) 



da^ 



da\f 



sr 



Q 



MC 

• 7(7 initial state 



aa ig 



( 



to -> tr 



(11) 



2m 2 



(12) 
(13) 



do i+ ' u) 



^ ^4 U) ^ 0) (4 U) ) <»*e ((4 U) ) 2 - e 



(u) 



.£(*) 



(14) 
(15) 
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• 7<Z initial state 



da^ q ,€} 



(0 



Y 1 %dzfP^{^) da qq Q ((,«)» - €«) , (16) 



a s d£, 



(a) 



2vr 



dz^P^(z^) da ig Q ((zi n) ) 2 - ei u) ) • (17) 



The Born cross sections da lgi da qq or deriyg have to be computed using the 
relevant definitions of s, t and u, as explained in ref. [9j. One also defines 13 1 



dsW - da ^ - ^ 7g 

™ " 19 u/t + t/U ~ 1 + P/u 2 



We stress that eqs. (|14p and (fT5|) describe QED branchings. When inserted in 
the MC@NLO short-distance cross sections, they subtract an MC contribution 
which is generated when the hadronic-photon component is showered. 

The situation of the contributions to the NLO expansion of PSMC cross 
sections is summarized in table [TJ 





qq -t QQ 


qq -t QQ 


19 QQ 


19 






-(t,u), Q(t,u), Q(t,u) 


10. 




+ («) 


-(*) 


iq 


+(*) 




-(«) 



Table 1: Short-distance contributions to MC subtraction terms, from Born processes qq — > 
QQ, qq — > QQ, and ■yg — ¥ QQ. The three rows correspond to real-emission processes, identified 
by their incoming partons. Each entry lists the emitting legs (+, — , Q, Q); for each emitting 
leg, we report in parentheses the different contributions I, according to the possible colour 
flows (which correspond to choosing = \T\/2). 



4. Comparisons to measurements of heavy- quark production at HERA 

In this section the MC@NLO predictions will be compared to selected mea- 
surements of various observables relevant to the production of B and of D*^ 
mesons, as reported by the HI and ZEUS experiments. The quark masses 
and PDFs used by MC@NLO are given in table [2] The uncertainty band 
in MC@NLO is computed by varying the factorization and renormalization 
scales independently, by a factor of two up and down from the default scale 
l/2(mq i T + ttlq t ), and by taking the envelope of the results obtained in this 
way. For charm quarks these scales may become fairly small, and indeed charm 
production is pushing the applicability of perturbative QCD to its limits - this 
is one of the reasons why it is interesting to compare MC@NLO results for this 
process to data. 
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4.75 GeV 


m c 


1.5 GeV 


Proton PDF 


Cteq6.6 


Photon PDF 


GRV 



Table 2: Parameter settings used by MCC8NLO. 



In order to quantify the level of agreement between MC@NLO predictions 
and data, the quantity ^ jndj has been calculated for each data set, by taking 
into account both theory and experimental uncertainties. The relevant results 
are summarized in table [S] for bottom measurements, and in table [5] for charm 
measurements . 

No effort has been made here to tune HERWIG (version 6.510) to HERA 
data, since the idea is that of making an out-of-the-box comparison. The only 
exception to this rule is an overall rescaling, applied to all hadron-level observ- 
ables in order to have values of branching ratios consistent with those reported 



by the PDG [14{ - these rescaling factors are equal to 1.34 and 1.5 for charm 



and bottom respectively. 

Jf.,1. B-hadron production 

Bottom-flavoured hadrons are typically tagged by searching for the muons 
that arise from If' s, which in turn come from the weak decay of the lowest-lying 
6-hadrons into lighter-quark states. These muons will in general have a large 
momentum transverse to the jet axis, so called p™ 1 . Also, the vertex from which 
these muons are radiated will be displaced relative to the hard interaction of 
the event, and this displacement is proportional to the lifetime of the 6-hadron. 
Often, only the transverse component 5 of this displacement is used in the b- 
tagging. These two methods of tagging the 6-quarks may be combined to further 
enhance the signal. 

In this section, comparisons will be made with three measurements per- 
formed by HI and ZEUS at HERA. They are: 

1. "Measurement of beauty photoproduction usin g de cays into muons in dijet 
events at HERA", by the ZEUS collaboration (lij; 

2. "Measurement of beauty production at HERA using events with muons 
and jets" , by the HI collaboration [16J ; 

3. "Beauty photoproduction measured using decays into muons in dijet events 



in ep collisions at *Js = 318 GeV", by the ZEUS collaboration (17 



These will be referred to as ZEUS-09, Hl-05 and ZEUS-03 respectively. The 
first two analyses use the combined method of both p™ 1 and 5 in the tagging of 
the &-quarks, while in ZEUS-03 only p™ 1 is used. The experimental cuts made 
for the bottom analyses are summarized in table [3] These cuts result in the 
visible cross-sections listed in table |4j where also the MC@NLO predictions are 
shown; theory and data are in good agreement. 



G 



Analysis 


ZEUS-09 


Hl-05 


ZEUS-03 




318 GcV 


318 GcV 


318 GeV 


Q 2 


< 1 GcV" 


< 1 GeV" 


< 1 GcV" 


vjb 


0.2 - 0.8 


0.2-0.8 


0.2-0.8 


Pt(M) 


> 2.5 GeV 


> 2.5 GeV 


> 2.5 GcV 




-1.6-1.3 


-0.55-1.1 


-1.48-2.3 


PT(jet 12 ) 


7, 6 GeV 


7, 6 GcV 


7, 6 GcV 


»7Cjet) 


-2.5-2.5 


-2.5-2.5 


-2.5-2.5 



Table 3: A summary of the cuts relevant to the bottom analyses considered in this paper. 



Visible a [pb] 


Measured 


MCQNLO 


ZEUS-09 


38.63:3 


42.08 ±4.91 


Hl-05 


38.4 ±6.38 


33.71 ±2.89 


ZEUS-03 


50.25 ±6.45 


48.39 ±3.87 



Table 4: The visible cross sections within the cuts listed in table[3] The different measurements 
are compared to MCCSNLO predictions. 



ZEUS-09 


MC@NLO 




0.18 


v{v) 


0.05 


Xj (jets) 


0.59 


A<Hjets) 


1.22 


A0(jets)|ir° bs < 0.75 


0.52 


Hl-05 


MC@NLO 




0.89 




0.11 


x 7 Qets) 


0.48 


ZEUS-03 


MCQNLO 




0.78 




0.34 


p T (b-jet) 


0.09 


p T (b) 


0.65 


Xj (jets) 


0.17 



Table 5: The x 2 / n< V f° r a U distributions in the bottom measurements shown in this paper. 



7 



> 

O 



> 

® -i 

o 1 



p, (b-Jet) [GeV] 
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Figure 1: Distributions of pT{b — jet) and Pt(V) from the measurement ZEUS-03. 
MC@NLO band includes full independent scale variations. 



The 



In fig. Q] the px of the jet containing a 6-quark is shown, as measured in 
ZEUS-03. This spectrum is then used by the ZEUS collaboration to reconstruct 
the pt spectrum of the 6-quarks, using the NLO calculation FMNR j3]. Both 
distributions are well described by MC@NLO. The pr{b — jet) prediction of 
MC@NLO, being at the hadron level, has been rescaled for the overall branch- 
ing ratio factor as discussed before. On the other hand, the pr(b) spectrum, 
being a quantity at the parton level, has not been rescaled. It should be stressed 
that the use of NLO computations matched to parton showers implies that the 
deconvolution of data from hadron to parton level is not necessary for a fair 
comparison with theoretical predictions; in fact, such deconvolutions have to 
be deprecated, since they introduce unnecessary theoretical biases in the mea- 
surements (e.g., in the present case, the underlying matrix elements of FMNR 
and MCQNLO are the same; hence, the comparison done for pr{b) is not as 
significant as that for pr(b — jet)). 

In fig. [5] the transverse momentum spectra of the tagged muons are shown. 
MC@NLO is describing all three data sets well for this observable. Also, the 
scale variations in MC@NLO are at the same level or smaller than the exper- 
imental uncertainties. The rapidity distributions of the muons are also well 
described by MC@NLO in all three data sets, as seen in fig. [3] 

Two types of correlations have been measured for bottom production. One 
is the x° hs distribution of the two leading jets in the measurements, shown in 
fig. HI The variable x° hs is defined by: 



a^Qeti.jeta) 



P_(jet)i+P_(jet 2 ) 



E 



All hadrons i 



P-ii) 



(19) 



where P_ = E - P z . Thus, x° bs is the fraction of the measured hadronic P_ 
carried by the two leading jets, and at the LO it coincides with the fraction 
of the photon energy which enter into the hard interaction. Therefore, the 
hadronic part of the calculation is expected to be more important for small 
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p.fri) [GeV] 
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Figure 2: Distributions of Pt(aO from the measurements ZEUS-09 (upper left), Hl-05 (upper 
right) and ZEUS-03 (bottom). The MC@NLO band includes full independent scale variations. 
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- ZEUS Data 

- MC@NLO 



10- 



-1.5 -1 -0.5 0.5 1 



n(M) 



j=r 



-0.4 -0.2 0.2 0.4 0.6 0.8 1 

1(H) 



— ZEUS DATA 
MC@NLO 



-1.5 -1 -0.5 0.5 1 1.5 2 

n(n) 



Figure 3: Distributions of r](fj,) from the measurements ZEUS-09 (upper left), Hl-05 (up- 
per right) and ZEUS-03 (bottom). The MC@NLO band is obtained with independent scale 
variations. 
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Figure 4: Distributions of x° bs from the measurements ZEUS-09 (upper left), Hl-05 (upper 
right) and ZEUS-03 (bottom). The MC@NLO band includes full independent scale variations. 

values of x° hs . MC@NLO describes all three measurements well. In these plots, 
the hadronic part of the MC@NLO calculation is also shown separate and, as 
expected, becomes significant for a;° bs < 0.75. 

An observable which is highly sensitive to higher-order effects is the differ- 
ence in azimuthal angle between the two leading jets. At Acf> ~ tt one observes 
the typical logarithmic divergence of infrared-sensitive variables computed at 
fixed-order in perturbation theory; this divergence is suppressed by the Sudakov 
damping, present in resummed computations and therefore naturally included in 
Monte Carlos. Multiple soft and collinear emissions are also relevant to small A<fi 
values (especially when rather small transverse momenta are involved), where 
however there are non-negligible contributions due to hard matrix elements, 
present in perturbative computations but not in ordinary PSMCs. It is there- 
fore clear that MC@NLO, by incorporating both contributions, is expected to 
give a better description for this variable than either fixed-order computations 
or ordinary PSMCs. We present the comparison of MC@NLO predictions with 
data in fig. [5j the agreement is satisfactory, given the large experimental uncer- 
tainties. This observable is also presented by separating the large and small s° bs 
regions, which are dominated by the pointlike and hadronic photon components 
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— ZEUS Data 
MC@NLO 



- ZEUS Data 

- MC@NLO 



A(Mjets) 



AiMjets) 



■ °!I0 2 - 



- ZEUS Data 

- MC@NLO 



A<Kjets) 



Figure 5: Distributions of A</>(jets), unbinned and binned in x° bs (jets) from ZEUS-09. The 
MC@NLO band includes full independent scale variations. 



respectively. The results are also displayed in fig. [5j 

4-2. D*^ production 

The D*^ mesons are detected through the so-called golden decay channel: 

D*± -> 2?°** K^^ ow . (20) 



The branching ratio for the golden decay channel is ~ 2.6% [14|, which is com- 
paratively low, but the advantage of this channel is that all the final state 
particles carry an electric charge, resulting in three charged tracks in the detec- 
tors. 

In this section MC@NLO will be compared to two D*^ measurements. These 
are: 

1. "Inclusive Z?*-Meson Cross Sections and D*-Jet Correlations in Photopro- 
duction at HERA" by the HI collaboration [18j; 

2. "Inclusive jet cross sections and dijet correlations in D*^ photoproduction 
at HERA" by the ZEUS collaboration 
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Analysis 


Hl-06 


ZEUS-05 




318 GeV 


318 GeV 


Q 2 


< 0.01 GeV 2 


< 1 GeV 2 


vjb 


0.29-0.65 


0.19-0.87 


Pt(D*) 


> 2 GeV 


> 3 GeV 


V (D*) 


-1.5-1.5 


-1.5-1.5 


?y(jet) 


-1.5-1.5 


-1.5-2.4 


pr(}et h2 ) 


4,3 GeV 





Table 6: A summary of the cuts in the D* meson measurements. 



Visible C-S [nb] 


Measured 


MC@NLO 


Hl-06 inclusive D* 


6.45 ±0.83 


6.45 ±0.78 


Hl-06 L>*+jcts 


3.01 ± 0.44 


2.88 ±0.29 


ZEUS-05 


6.80 ±0.26 


5.77 ±0.42 



Table 7: The resulting visible cross-sections from the cuts listed in table [6] for the different 
measurements as well for the MC@NLO predictions. Hl-06 consists of two measurements: of 
inclusive D* and of D* + jets. These are listed separately here. 

These will be referred to as Hl-06 and ZEUS-05 respectively. The experimental 
cuts applied in the D* analyses are summarized in table [6l These cuts result 
in the visible cross sections reported in table [3 together with the theoretical 
predictions. MC@NLO describe the Hl-06 measurements very well and is two 
sigma below the data for ZEUS-05. 

4-2.1. Inclusive D*^ production 

In fig. [6] the px(D*) and rj(D*) spectra are shown for inclusive D* measure- 
ment. The Hl-06 measurement is well described but the scale uncertainties in 
the MC@NLO prediction are rather large in comparison with the experimental 
uncertainties. In the pr(D*) distributions, these scale uncertainties are reduced 
for larger pt(D*), since the overall hardness of the process is increased. In fig. [5] 
the r](D*) distribution is also shown separately for pt{D*) > 4.5 GeV, where 
the scales uncertainties in MC@NLO are at the same level as the experimental 
uncertainties. Here, two of the bins are well described, but the overall shape is 
not. However, given the very large uncertainties involved, the x 2 / n( V for this 
distribution is a mere 1.67. 

4-2.2. D*^ plus jet production 

When demanding the tagging of a hard jet, as well as that of a D*^ meson, 
the scale dependencies in MC@NLO are expected to be smaller than for inclusive 
variables, owing to the extra hardness of the jet. In fig. [3 the pr(D*) and rj(D*) 
distributions from the Hl-06 measurements are shown, when events with a hard 
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- HI Data 

- MC@NLO 



-1.5 -1 -0.5 



ti(D*) 



- HI Data 

- MC@NLO 



4.5<p (D*)[GeV]<8.5 



-1.5 -1 -0.5 



1 1.5 
T1(D*) 



Figure 6: Distributions of px{D*) and rj(D*) of inclusive D* measurements Hl-06. The 
MC@NLO band includes full independent scale variations. 
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Hl-06 inclusive 


MC@NLO x 2 /nds 


Pt(D*) 


1.41 


V(D*) 


0.46 


r](D*) p T (D*) > 4.5 GeV 


1.67 


Hl-06 L>*+jct 


MC@NLO x 2 /nds 


Pt(D*) 


2.49 


T)(D*) 


0.52 


obs 
•by 


1.94 


A<£(£>*,jet) 


0.67 


A V (D*,)et) 


0.28 


ZEUS-05 


MC@NLO x 2 /nds 


Pt(D*) 


0.94 


j,ohs 


1.05 


A<^D*,jet) 


1.37 


PT(jj) 


1.01 


Ma 


1.42 


^(untagged jet) pr{jet) > 9 GeV 


0.78 



Table 8: The x 2 / n df for all distributions in D'^-measurements shown. 



jet with pr(jet) > 4 GeV are chosen. The scale dependence is indeed reduced. 
The data are still well described, even though the pr{D*) spectrum is a bit 
harder in MC@NLO than in the Hl-06 data, while in ZEUS-05 Q the opposite 
is observed. 

Also, correlations between D* mesons and jets have been measured in Hl- 
06 and ZEUS-05. In fig. distributions in x° hs are shown. One can see that 
the contribution by the hadronic part of MC@NLO is larger in the Hl-06 mea- 
surement. However, MC@NLO is one sigma above the data for large x° hs in 
this measurement, while the whole spectrum is well described in the ZEUS-05 
measurement. 

In fig. IH1 the difference in azimuthal angle between the D* meson and the 
hardest jet not containing the D* (for Hl-06), and between the two hardest jets 
in a D* event (for ZEUS-05) are shown. One can see that MC@NLO describes 
the data over the whole A<f) spectrum for both analyses, something fixed order 
NLO calculations cannot do (see e.g. fig. 9 in ref. [3], and fig. 11 in ref. (l9l|). 

Another observable which is sensitive to higher order effects is the pr of the 
two leading jets, which was measured in ZEUS-05, and shown here in fig. 1101 In 
this figure, the distribution of the invariant mass of the two jets is also presented. 



1 This pt(D*) distribution is not presented in ref. [lSj - We have compiled it from other 
distributions binned in prp{D*). 
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> 

O 



- HI Data 

- MC@NLO 



~ 1f 
Q 



10'^ 



9 10 11 

p, (D*) [GeV] 



10 12 14 16 18 20 

p(D*) [GeV] 



ti(D*) 



Figure 7: Distributions of pt(D*) and r](D*) in the D* + jets measurements from Hl-06 as 
well as pt(D* ) from ZEUS-05. The MCONLO band includes full independent scale variations. 




Figure 8: Distributions of x° bs (jets) from Hl-06 (left) and ZEUS-05 (right). The MCONLO 
band includes full independent scale variations. 
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20 40 60 80 100 120 140 160 180 0.5 1 1.5 2 2.5 3 

A<KD*, jet) p] A«MIJ) 

Figure 9: Distributions of A0 from Hl-06 (left) and ZEUS-05 (right). The MCONLO band 
includes full independent scale variations. 

MC@NLO is seen to describe both these observables in a reasonable manner. 
Also, the difference in pseudo-rapidity between the D* and the hardest other jet 
is sensitive to higher order radiations. This was measured in Hl-06 as shown in 
the bottom plot of fig. [TO] MC@NLO is describing the spectrum in Ary(_D*,jet) 
very well, with x 2 /W/ = 0.28. 

5. Conclusions 

In this letter, we have applied the MC@NLO formalism to the pointlike 
photoproduction of heavy-quark pairs. Together with the already- available QQ 
hadroproduction code, this has been employed to carry out a comparison be- 
tween theoretical predictions and the data collected by the HI and ZEUS ex- 
periments at HERA, relevant to b- and c-hadron observables. 

In particular, MC@NLO has been compared to five sets of measurements, 
three for ^-flavoured hadrons, and two for c-fiavoured hadrons. All data have 
been shown to be described within one standard deviation by MC@NLO. It 
should be pointed out that MC@NLO predictions are absolute, except for an 
overall rescaling needed to obtain a branching ratio into muons (in the case of 
bottom production), and into Z?*'s (in the case of charm production) compatible 
with the world averages. 

Although the overall agreement between theory and data appears to be satis- 
factory, it should be kept in mind that the uncertainties involved are sometimes 
quite large. For bottom production, theoretical errors are at the same level 
or smaller than experimental ones. On the other hand, for charm production 
the largest uncertainties are those of MC@NLO, and therefore the comparisons 
carried out here only loosely constrain perturbative QCD predictions. The sit- 
uation improves when hard jets are also part of the observable definitions. 
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Figure 10: Distributions of p T (jj) and Mjj from ZEUS-05 and of Ar;(D*, jet) from Hl-06. The 
MC@NLO band includes full independent scale variations. 
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